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During the first 500 million years of cosmic history, the first stars 

and galaxies formed, seeding the Universe with heavy elements and 
eventually reionizing the intergalactic medium [1-3]. Observations 
with JWST have uncovered a surprisingly high abundance of can- 
didates for early star-forming galaxies, with distances (redshifts, 
z), estimated from multi-band photometry, as large as z = 16, 
far beyond pre-JWST limits [4-9]. While generally robust, such 
photometric redshifts can suffer from degeneracies and occasion- 
ally catastrophic errors. Spectroscopic measurement is required to 
validate these sources and to reliably quantify physical properties 
that can constrain galaxy formation models and cosmology [10]. 
Here we present JWST spectroscopy that confirms redshifts for two 
very luminous galaxies with z > 11, but also demonstrates that 
another candidate with suggested z = 16 instead has z = 4.9, with 
an unusual combination of nebular line emission and dust redden- 
ing that mimics the colors expected for much more distant objects. 
These results reinforce evidence for the early, rapid formation of 
remarkably luminous galaxies, while also highlighting the necessity 
of spectroscopic verification. The large abundance of bright, early 
galaxies may indicate shortcomings in current galaxy formation 
models, or deviation from physical properties (such as the stellar 
initial mass function) that are generally believed to hold at later 
times. 

The spectroscopic targets were selected from JWST/NIRCam [11] imag- 
ing obtained through the Cosmic Evolution Early Release Science (CEERS) 
Survey [12] because their colors at wavelengths of 1 — 5 microns matched 
predictions for galaxies with redshifts z > 9 [9], including the key signature 
of a “break” in flux density at wavelengths shorter than that of redshifted 
hydrogen Lyman a (Lya; rest-frame wavelength Arest = 1216 A) due to the 
largely-neutral intergalactic medium (IGM) expected at early times. Here we 
focus attention on three bright extreme-redshift candidates: two with photo- 
metric redshifts z ~ 11 (CEERS2_5429, aka Maisie's Galaxy [5, 9], here MSA 
ID 1; and CEERS2_588 [9], here MSA ID 10), and another with a photomet- 
ric redshift z ~ 16.4 (CEERS-93316 [6]; here MSA ID 0). At their estimated 
redshifts, these three galaxies would have UV luminosities ~ 2 to ~ 9 times 
brighter than any other spectroscopically confirmed galaxy at z > 11 [13, 14]). 
We show these photometric redshift estimates in Figure 1. 

Although independent analyses by different teams [5-7, 9, 15] identified 
these candidates using 7-band JWST/NIRCam photometry, it is possible for 
foreground objects to infiltrate high-redshift samples due to unusual colors 
resulting from older stellar populations, strong nebular emission lines, and/or 
dust reddening [6, 16-18]. A tentative detection at 850 um of dust emis- 
sion in the vicinity of CEERS-93316 [18] could support the dust reddening 
hypothesis, but we have investigated this using new observations at 1.1 mm 
with the NOEMA interferometer and localized the dust emission to another 
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Figure 1: Photometric redshift probability density function for the 
primary high-z candidates. The likelihood (bottom) and x? goodness-of- 
fit (top) distribution functions for the photometric redshift, derived from the 
JWST CEERS photometry. Results for Maisie's Galaxy, CEERS2.588 and 
the z ~ 16.5 candidate galaxy CEERS-93316 are shown in blue, green and 
red, respectively. The line shading denotes the results using the photometric 
measurements from each team. The spectroscopic measurements of z = 11.416, 
z = 11.043 and z = 4.912 from this paper are indicated by the thick vertical 
lines. 


galaxy approximately 1.5 arcsec away (see Methods for details). Nonetheless, 
reddening could still affect the observed colors of CEERS-93316. 

We observed these galaxies with JWST/NIRSpec [19] using its multi- 
object spectroscopy mode and prism disperser, whose low spectral resolution 
at wavelengths shorter than 2.5 jum affords high sensitivity for detecting faint 
continuum emission and spectral breaks. Details of the observations and data 
analysis, as well as results from simultaneous observations of several other 
candidate z > 9 galaxies, are summarized in the Methods section. 

The spectrum of Maisie's Galaxy (Figure 2) shows well-detected contin- 
uum emission with a sharp break, below which no significant flux is detected. 
The most likely interpretations of this feature are the Lya IGM break or the 
Balmer break (Arest = 3646 A); the latter is prominent in galaxies where star 
formation ceased more than 100 Myr before the time of observation. Integrat- 
ing the spectrum in wavelength intervals redder and bluer than the observed 
break, we measure a 3c limit on the flux density ratio > 2.6, inconsistent 
with values expected for the Balmer break [14] (see Methods). We used several 
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Figure 2: Color stamps and spectra of the confirmed z > 11 galaxies. 
Panels a, c) Color images of Maisie’s Galaxy (MSA ID 1) and CEERS2_588 
(MSA ID 10), showing approximately the true rest-frame ultraviolet color 
(made from the F277W, F356W, and F444W images). Panels b, d) Spectra of 
Maisie’s Galaxy and CEERS2_588 highlighting the detected Lya breaks and 
[O 11] emission lines. The top panel shows the 2D spectrum. The red horizontal 
lines show the spatial extent over which the 1D spectrum, shown in the bottom 
panel, was extracted. The shaded grey area corresponds to 1o uncertainties in 
the 1D spectrum and the zero flux level is marked by a dashed black horizon- 
tal line. A clear spectral break is seen, below which there is no significant flux. 
We interpret this as the Lya break at a redshift of 11.416 and 11.043, which 
corresponds to a time 392 and 410 Myr after the Big Bang, respectively. 


methods [20] to measure the wavelength of the spectral break (see Meth- 
ods and Extended Data Table 1). Interpreting this feature as the Lya break 
(Figure 2b) yields a redshift z = 11.44*0 04, consistent with photometric red- 
shift estimates. We also identify an emission line at 4.629 um as the unresolved 
singly-ionized [O 11] doublet (Arest = 3727,3730 A) at z = 11.416. 

CEERS2.588 exhibits a comparably strong Lya break and [O 11] emission 
line at z = 11.043 [21]. Its break amplitude is > 2.2, again inconsistent with 
a Balmer break. Another, fainter galaxy (MSA ID 64) also has a strong (> 
1.9) break and no emission lines, consistent with a Lya break redshift z = 
10.10*0:3 (see Methods). 

NIRSpec observations of star-forming galaxies at z « 9.6 frequently detect 
strong [O ri] (Arest = 4960,5008 A) and HB (Arest = 4863 A) emission lines 
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Table 1: Summary of galaxy properties 


Property MSA ID 1 MSA ID 10 MSA ID 0 
(Maisie's Galaxy) (CEERS2.588)  (CEERS-93316) 
RA (J2000) 214.943152 214.906640 214.914550 
Dec (J2000) 52.942442 52.945504 52.943023 
Photometric Redshift 11.08+9:36 Hour d 16,45" 0: 
Spectroscopic Redshift 11.416+0.005 11.043+0.003 4.912+0.001 
TpigBang (Myr) 392 410 1194 
Muv -201491 -—20.3+2:1 —16.2*9.5 
log( M. / Mo) 8.6 E 0.3 8.7 E 0.1 8.80.1 
Range: 8.0 — 9.0 8.1 — 9.1 8.3 —9.2 
Ay (mag) 0.1+0.1 0.2+0.1 2.3+0.2 
Range: 0.0 — 0.8 0.03 — 0.7 0.8 — 2.7 
SFR (Mo yr7') 241 104 60 + 20 
Range: 1-12 2-17 2 — 80 
log(sSFR/yr-!) -83*02 2 ari emp 
Range: -8.6 — -7.6 -8.2 — -7.5 -8.1 — -6.9 


Table 1: The photometric redshift calculation is described in the Methods 
and has been updated for this work. Spectroscopic redshifts are derived from 
emission lines; see Methods for a full discussion. Tpigpang denotes the age of 
the universe at the spectroscopic redshift. All uncertainties are 68% confidence 
intervals. Stellar masses (M,), SFRs averaged over the last 100 Myr and Ay 
values have been derived using three different stellar population modeling tech- 
niques (see Methods). For each property, the first line gives the best-fit values 
and 68% confidence intervals from the method simultaneously fitting both the 
spectra and the photometry (Synthesizer), while the second line shows the 
ranges spanned by all three modeling methods including their 68% confidence 
intervals. 


[22-24] that shift beyond the instrument's spectral coverage for redshifts z > 
9.6. Emission lines at blue optical and ultraviolet rest-frame wavelengths are 
typically weaker and harder to detect in low signal-to-noise ratio (S/N) spectra. 
The [O ri] lines for Maisie's Galaxy and CEERS2_588 are only weakly detected 
(4.90 and 6.40, respectively) even with moderately long NIRSpec exposure 
times, highlighting the importance of a robust Lya break detection. 

Maisie's Galaxy and CEERS2.588 are the most UV-luminous galaxies with 
confirmed spectroscopic redshifts z > 11 to date, with rest-frame UV abso- 
lute magnitudes Myy = —20.1 and -20.3, respectively. The physical properties 
derived through stellar population modeling (see Methods for details) on the 
observed photometry and spectroscopy of Maisie's Galaxy and CEERS2.588, 
fixing the redshifts to the spectroscopic values, are presented in Table 1. We 
find that they have stellar masses log;4(M,/Mc) = 8.6 to 8.7 and low dust 
extinction (Ay = 0.1 to 0.2 mag), with high star-formation rates (SFR) for 
their stellar masses (log; (sSFR = SFR/M,) = —8.3 to -7.7). These phys- 
ical properties are consistent with similar values found for the few other 
spectroscopically confirmed galaxies at z > 11 [13, 25]. 
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Figure 3: NIRSpec spectrum of Zphot ~ 16 candidate CEERS-93316 
(MSA ID 0). Panel a) shows the 2D spectrum with the red lines denoting 
the spatial extraction limits, and panel c) shows the 1D extraction. Emis- 
sion lines labeled in panel c) indicate an unambiguous z = 4.912 for this 
object. The shaded grey area corresponds to lo uncertainties in the 1D spec- 
trum. Panel c) also illustrates the wavelengths covered by the JWST/NIRCam 
filter set used to select this galaxy, highlighting how [Om] contamination 
in F277W (green) and Ha contamination in all of F356W (yellow), F410M 
(orange) and F444W (red) conspired to mimic the signature of both a very 
high-redshift Lyo break and a blue continuum at longer wavelengths. This 
unusual (but educational) situation is detailed further in Figure 4. Panel b) 
shows a F150W+F200W+F277W color image of this galaxy. 


The spectrum of CEERS-93316 (Figure 3) shows emission lines from the 
hydrogen Balmer series, [Om], [Om], and [St] (Arest = 6718,6733 A) at a 
redshift z = 4.912 + 0.001. This is the same redshift as for another galaxy, 
CEERS-DSFG-1 (MSA ID 2) [18, 26], which lies 26 arcsec away and was pre- 
viously discussed as a dusty, star forming galaxy with colors that could mimic 
a high-redshift Lya break [18]. At least two other neighboring galaxies have 
very similar redshifts (see Methods). The continuum emission of CEERS-93316 
fades at wavelengths shorter than 2.5 jum (rest frame SS 0.4 um), although flux 
is significantly detected in the spectrum down to A z 1.7 jum, and marginally 
detected down to A ~ 1.3 um. The strong Ha+[NI]] and [O r11]--H lines have 
large rest-frame equivalent widths (EWs; 760 + 67 A and 1116 + 13 A, respec- 
tively), and hence contribute significantly to the flux measured by NIRCam 
through the F277W, F356W, F410M, and F444W filters. At this precise red- 
shift, Ha emission actually contributes to the flux in the last three of those 
filters (Figure 3). The combination of the strong lines and the red continuum 
leads to the red colors measured between filters F150W, F200W, and F277W 
and the blue colors from F277W to F444W (Figure 4c). 

The significant reddening implied by the faint, red UV continuuum emission 
is corroborated by the large Balmer decrement measured from the spectrum 
(Ha/H8 = 5.2 + 0.5), which implies — despite the large EWs of the lines — 
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Figure 4: Spectral energy distributions of the three main sources 
reported. For a) Maisie’s Galaxy (MSA ID 1), b) CEERS2_588 (MSA ID 10) 
and c) CEERS-93316 (MSA ID 0), each panel shows the observed photometry 
(colored points) compared to the best-fitting stellar population model (black 
line) and the observed spectra (colored line). In panel c), the model is delib- 
erately shown at z = 16.2 to demonstrate the way in which strong emission 
lines impact the observed photometry in this specific case to replicate the pres- 
ence of the Lya break at z ~ 16 and an apparently blue continuum at longer 
wavelengths. The brown squares show the observed photometry omitting the 
spectroscopic emission line contribution (shaded brown in the spectrum). Error 
bars are lo uncertainties while the upper limits correspond to 2c limits. 
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a substantial nebular color excess of E(B — V) = 0.60 + 0.10 (see Methods). 
Indeed this combination of line emission and reddening at z ~ 5 has been con- 
sidered previously as an alternative explanation for the colors of CEERS-93316 
[6, 16-18]. The Ha emission may imply rapid star formation (28*5 Mo yr-!, 
see Methods), although some other faint red objects at z ~ 5 show strong 
line emission with broad Balmer line components that indicate the presence 
of active galactic nucleii [e.g., 27, 28]. 

The example of CEERS-93316 is a reminder that the prevalence of very 
strong nebular line emission in early galaxies can influence broad-band pho- 
tometry with instruments such as JWST/NIRCam and mimic the colors 
predicted for much more distant objects. Exceptional candidates for galaxies 
at extreme redshifts require spectroscopic confirmation. However, we note that 
this is a fairly extreme case, as it is only at very specific redshifts that strong 
emission lines can fully mimic a Lyo break with multiple broad-band filter 
detections, and reddening is also required to suppress detection in bluer filters 
(see Methods). 

Fortunately JWST/NIRSpec is a powerful tool for spectroscopic analysis 
and redshift measurement, as demonstrated by our verification of two galaxies 
with z > 11 and three other galaxies at 7.9 < z < 10.1 (discussed in Methods). 
NIRSpec observations have now confirmed redshifts for 10 galaxies at z > 10 
[14, 20, 21, 29, 30], including five with Muy < —20, unexpectedly bright 
compared to most pre-J WST predictions. These confirmations are crucial to 
validate the methods employed for z 2 10 photometric selection [20] and to 
improve our physical understanding of the first phases of galaxy formation. 

Early analyses of candidate galaxies at z = 9.5-12 have shown that their 
observed abundance is significantly higher than predicted by nearly all theoret- 
ical models [5, 7, 9, 31]. This could indicate that specific assumptions in these 
models need to be modified at very high redshifts, including the possibility 
of negligible dust attenuation [32, 33] or an increased efficiency of converting 
gas into stars, or even a shift toward the formation of more massive stars in 
these early epochs. Such a “top-heavy” initial mass function (IMF) has indeed 
been predicted to be present in early galaxies, driven by a combination of 
their chemically simple makeup [34, 35] as well as a rising cosmic microwave 
background temperature floor [36]. However, before any such unusual phys- 
ical properties can be robustly inferred for early galaxies, their photometric 
redshift estimates must be verified. Four out of seven photometric candidates 
from CEERS data at 9.5 < z < 12 with Myy < —20 have been spectro- 
scopically confirmed [20], including the present examples at z > 11, with no 
confirmed interlopers except the z ~ 16 candidate CEERS-93316. The space 
density of confirmed objects alone exceeds the predictions of most theoretical 
models. The spectroscopic observations reported here reinforce the higher- 
than-expected abundance of such bright galaxies in the very early universe, 
while also vetting photometric redshift outliers like CEERS-93316 that could 
be misinterpreted if left unrecognized. 
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Methods 


Cosmological model 


Throughout this work we adopt a A cold dark matter cosmological model with 
parameters Qm = 0.315, Q4 = 0.685 and Ho = 67.36 km s^! Mpc^! [37]. All 
magnitudes are presented in the AB system [38]. 


Near-infrared imaging data and spectroscopic target 
selection 


We select high redshift galaxy candidates from the first epoch of imaging 
from the Cosmic Evolution Early Release Science (CEERS; PI: S. Finkelstein 
[5]) survey. These data consist of four roughly contiguous pointings with the 
NIRCam instrument on JWST, covering an area of ~35 arcmin?. The full area 
is covered by six broad-band filters (F115W, F150W, F200W, F277W, F356W, 
F444W) and one medium-band filter (F410M). We use the v0.5 data release 
from the CEERS team! [12]. These reductions use the STScI JWST pipeline 
with several custom modifications, correcting for various issues, including 1/ f 
noise, wisps, and snowballs?, with updated photometric calibration from late 
2022. The 50 photometric depths of these data for compact source detection 
are ~29-29.2 AB magnitudes in F115W, F150W, F200W, F277W and F356W, 
and —0.7 and —0.5 mag shallower in FA10M and F444W, respectively. 

Photometry was optimized for accurate measures of colors, generally fol- 
lowing previous CEERS measurements [9], with a few key updates to improve 
measured colors, particularly between HST and JWST. Specifically, all bands 
with point spread functions (PSFs) smaller than that of F277W (i.e., ACS 
F606W, F814W, and NIRCam F115W, F150W, and F200W) were convolved to 
match the F277W PSF. For the remaining images, fluxes were measured in the 
native images, but a correction was applied on a per-source basis as the ratio 
of the flux in the native F277W image to that of the F277W image convolved 
to match the PSF in a given larger-PSF image. We then used simulation- 
based corrections to derive accurate total fluxes. Photometric redshifts were 
derived using EAZY [39], including new custom templates designed to match 
the expected blue colors of early galaxies [40]. 

The NIRSpec spectroscopic observations targeted three primary galaxy 
candidates: CEERS2.5429, aka “Maisie’s Galaxy" [5, 9], with photometric 
redshift z ~ 11.1, “CEERS-93316” with photometric redshift z ~ 16.4 [6], 
and “CEERS-DSFG-1”, a JWST-detected galaxy with millimeter-wavelength 
emission that is likely to arise from dust-obscured star formation at z ~ 5 
[26], but whose near-infrared colors could be interpreted as suggesting z > 10 
[18]. Many other faint, high-redshift candidates were also observed simultane- 
ously, including CEERS2_588, which has photometric redshift z ~ 11 [6, 9] and 
is comparably bright as Maisie's Galaxy. These object names come from the 


lhttps:/ /ceers.github.io 
?https://jwst-docs.stsci.edu/data-artifacts-and-features/snowballs-and-shower-artifacts 
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original papers that reported these objects rather than from a homogeneous 
source catalog. 


NOEMA data and analysis 


Observations with the JCMT and its SCUBA2 bolometer array detected ten- 
tative 850 um emission in the vicinity of CEERS-93316. Previous analysis [18] 
showed that this submillimeter emission is unlikely to arise from an object at 
z > 10, and that a z ~ 5 galaxy with dust attenuation and strong emission 
lines could explain the measured JWST photometry. 

We observed dust continuum with band 3 of the Institute Radio Astronomie 
Millimétrique (IRAM) NOrthern Extended Millimeter Array (NOEMA) as a 
Director's Discretionary Time (DDT) program #D22AC001 (PI: S. Fujimoto). 
The observations were carried out between 2022 Nov 26 and 2022 Dec 08 in 
several visits with the C configuration using 10-11 antennas. The data were 
processed in the standard manner with the pipeline using the latest version of 
GILDAS. We used CASA version 6.5 for the imaging [41]. 

The upper and lower side bands of the band 3 receiver were tuned to 
256.6 GHz and 271.1 GHz, resulting in a central wavelength of 1.134 mm. 
The observations of 3C45, 3C273, J1439+499, and 3C84 served as band-pass 
phase calibrators. Additional targets 1418+546 and J1439+499 were used for 
the phase and amplitude calibrations. We calibrated the absolute flux scale 
against MWC349, LkHal01, and 2010+723 whose flux densities are regularly 
monitored at NOEMA. The total integration time on-source was 6.5 hours. 

To maximize sensitivity, we used natural weighting for the imaging. The 
resulting 1 mm continuum map has synthesized beam FWHM of 0.791 x 0."67 
with lo sensitivities of 28 Jy beam! for the continuum. We do not identify 
any line features at the target position and thus produce the 1 mm continuum 
map from all channels. 

In Extended Data Figure 1, we show the NOEMA 1.1-mm intensity con- 
tours overlaid on the NIRCam color image around CEERS-93316. We do not 
detect emission at the position of CEERS-93316, with a 3e upper limit of 
0.084 mJy beam !. 

On the other hand, the NOEMA observations revealed relatively bright 
continuum emission at 9c significance from a nearby (~ 1.5") galaxy with 
a photometric redshift z ~ 4.9, very similar to the spectroscopic redshifts 
of CEERS-93316, CEERS-DSFG-1 (also detected at millimeter wavelengths), 
and other several galaxies in the field (see below and Extended Data Table 1). 
We measure a 1.1 mm flux density 0.56 + 0.07 mJy with an optimized aper- 
ture. This is in good agreement with the previous report of SCUBA2 850 um 
emission with 0.64 + 0.26 mJy [18], implying that this neighboring galaxy is 
the main contributor to the SCUBA2 detection. 
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A 


CEERS-93316 


Extended Data Figure 1: NOEMA 1.1-mm contours overlaid on the NIR- 
Cam 6" x 6" RGB color image (R: FA44W, G: F356W, B: F200W) around 
CEERS-93316 (MSA ID 0) pointed by the white arrow. North is up and East 
is left. The magenta contours are drawn at lo intervals from +20 to +9o. 
CEERS-93316 is undetected at 1.1-mm, while robust emission from a nearby 
source is likely responsible for a previously reported tentative SCUBA2 detec- 
tion [18]). 


Near-infrared spectroscopy and data reduction 


JWST/NIRSpec multi-object spectroscopic observations (DDT program 
#2750; PI: P. Arrabal Haro) were executed on UTC 2023 March 24-25 using 
the PRISM/CLEAR combination of disperser and filter. The prism disperser 
covers the wavelength range 0.6 to 5.3 um with varying spectral resolution 
R= A/AA & 30 at A = 1.2 um to > 300 at A > 5m. The low resolution of 
the prism at bluer wavelengths aids the detection of faint UV continuum and 
the Lya break, while the higher resolution at red wavelengths facilitates the 
detection of redshifted optical rest-frame emission lines. 

The NIRSpec Micro-Shutter Assembly (MSA) [42] was configured to 
observe 160 faint-object targets. MSA ID numbers were assigned to the tar- 
gets when the spectroscopic observing plan was designed. For most observed 
targets, three-shutter slitlets were employed, enabling a three-point nodding 
pattern to facilitate background subtraction. Exceptionally, CEERS2.588 was 
fit into the MSA configuration *by hand" using a two-shutter slitlet, and there- 
fore received two thirds of the full exposure time. The observations consisted 
of 3 sequences of 3 nodded integrations, each with 9 to 10 readout groups in 
NRSIRS2 readout mode, for a total combined exposure time of 18387 s. 

The MSA configurations included a total of 7 NIRCam-selected z = 8.5 
galaxy candidates as well as 153 other targets. We report results for 10 objects 
in Extended Data Table 1 including the primary targets described above. 
Essential parameters for Maisie's Galaxy (MSA ID 1), CEERS2.588 (MSA ID 
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Extended Data Table 1: Summary of spectroscopically confirmed galaxies 
MSA R.A. Dec F277W phot pings abreak zeont 
ID (J2000) (J2000) (AB) 68% C.L. 

1 214.943152 52.942442 27.9 10.7-11.5 11.416 + 0.005 11.44100 11.3870 18 
10 214.906640 52.945504 27.3  10.6-11.3 11.043 + 0.003 11.4343 11.031902 
64 214.922787 52.011529 28.3  10.4-11.6 — 10.10*013 9.737021 


28 214.938642 52.911749 26.9 8.90-9.00 8.763 + 0.001 8.907005  8.72*0 01 
355  214.944766 52.031450 28.7 8.0-8.5 7.925 + 0.001 7.717042 7.8875 01 
0 214.914550 52.943023 26.5  16.0-16.6 4.912 + 0.001 — 4.89* 0-01 
2 214.909116 52.937205 26.2  3.4-13.6 4.910 4 0.003 — 4.88* 0-01 
2763 214.927789 52.935859 25.1 4.87-5.44 4.902 + 0.001 = LRITUOE 
3149 214.914917 52.943622 24.7 0.88-4.81 4.901 + 0.001 = ABT? 01 


69 214.861602 52.904604 28.2 9.5-10.3 — = E 
Extended Data Table 1: The photometric redshifts limits correspond to 68% 
confidence intervals. Spectroscopic redshifts presented are based on fitting of 
emission lines, the Lya break, and the complete continuum measured with the 
NIRSpec prism. Maisie's Galaxy, CEERS2.588 and CEERS-93316 correspond 
to MSA IDs 1, 10 and 0, respectively. 


10), and CEERS-93316 (MSA ID 0), including results from the present anal- 
ysis, are summarized in Table 1. Detailed analysis of CEERS-DSFG-1 (MSA 
ID 2) will be presented elsewhere. 

'The NIRSpec data processing followed the methodology employed for other 
CEERS NIRSpec prism observations [20, 22]. We make use of the STScI Cal- 
ibration Pipeline? version 1.8.5 and the Calibration Reference Data System 
(CRDS) mapping 1088. We use the calwebb_detector1 pipeline module to 
subtract the bias and the dark current, correct the 1/f noise, and generate 
count-rate maps (CRMs) from the uncalibrated images. The parameters of 
the jump step are modified for an improved correction of “snowball” events 
associated with high-energy cosmic rays. 

The resulting CRMs are then processed with the calwebb spec2 pipeline 
module which creates two-dimensional (2D) cutouts of the slitlets, performs the 
background subtraction making use of the nodding pattern, corrects the flat- 
fields, implements the wavelength and photometric calibrations and resamples 
the 2D spectra to correct the distortion of the spectral trace. The pathloss 
step accounts for the slit loss correction at this stage of the reduction process. 

The images of the individual nods are combined at the calwebb spec3 
pipeline stage, making use of customized apertures for the extraction of the 
one-dimensional (1D) spectrum. 

The 2D and 1D spectra are simultaneously inspected with the Mosviz visu- 
alization tool* [43] to mask possible remaining hot pixels and other artifacts 
in the images, as well as the detector gap (when present). Data from the three 
consecutive exposure sequences are combined after masking image artifacts to 
produce the final 2D and 1D spectral products. 


3https://jwst-pipeline.readthedocs.io/en/latest/index.html 
^https://jdaviz.readthedocs.io/en/latest / mosviz/index.html 
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The JWST pipeline uses an instrumental noise model to calculate flux 
errors for the extracted spectra. We test and rescale these flux errors for 
the effect of data correlation introduced by the pipeline following procedures 
described elsewhere [20]. To further test possible flux discrepancies due to 
uncorrected slit losses or inaccurate flux calibration, we integrate the spectra 
through the NIRCam filter bandpasses and compare the resulting synthetic 
photometry to the measured NIRCam photometry. For the main sources pre- 
sented in this work, the flux ratios are similar for all NIRCam bands with 
robust (5c) detections, and do not show any trend with wavelength. We rescale 
the spectroscopic fluxes by a constant factor derived from the average of the 
flux ratios for all bands with detections. These rescaling factors are 1.29, 1.46 
and 1.56 for CEERS-93316 (MSA ID 0), Maisie's Galaxy (MSA ID 1) and 
CEERS2_588 (MSA ID 10), respectively. 


Redshift estimation 


For objects with multiple clear emission lines, we were easily able to identify the 
[O 11] and/or Ha emission lines. Spectroscopic redshifts are measured by com- 
paring the observed emission lines with their theoretical vacuum transitions. A 
Gaussian profile is fitted for every spectral line to establish the transition wave- 
length. The Gaussian fluxes are used to weigh the mean redshifts in Extended 
Data Table 1. In those observations, where the [O 111] doublet lines are merged, 
the profile fitting assumes the same gas kinematics and their amplitudes fixed 
by the theoretical emissivity ratio (2.98 [44]) to isolate them. 

We also estimate the redshifts using the detected continuum with two meth- 
ods. The first method consists of fitting a five-parameter model to describe an 
object’s spectrum and determine the wavelength of the Lya break [20]. We 
first create the intrinsic spectrum based on the redshift, UV absolute mag- 
nitude, and UV spectral slope. We then add IGM absorption by zeroing out 
the spectrum below 1215.67 A rest-frame wavelength, and add Lya damping 
wing absorption by adding two additional parameters: the neutral hydrogen 
fraction (ayy) and an ionized bubble radius (Rpubpie) [14, 45]. We derive pos- 
terior constraints on these five parameters using an IDL implementation of the 
EMCEE Python code [46], restricting the spectrum to wavelengths below 2500 
A rest-frame for a given redshift. The second method is full spectral fitting 
to the NIRSpec data, which detects the continuum for all our galaxies, using 
EAZY [39] and a variety of spectral templates, including UV-bright and dusty 
galaxies. For Maisie’s Galaxy and CEERS2_588 the spectra cover both the 
Lyman and Balmer breaks. All methods provide results in agreement within 
Az ~ 0.1, consistent within our estimated uncertainties. 


Redshift results 


We derive redshifts from multiple emission lines for galaxies with MSA IDs 
0 (via [On], [Om], Ha), 2 (via [Orr], Ho), 28 (via [On], [Om]) and 355 
(via [O rt], H8). We measure Lya break redshifts for MSA IDs 1, 10, 28, 64 
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Extended Data Figure 2: 2D and 1D NIRSpec spectra of six high-z candi- 
dates. Features and spectroscopic redshifts are indicated, except for MSA ID 
69 where we cannot determine a redshift with confidence. The shaded grey 
area corresponds to 1c uncertainties. 
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(Maisie's Galaxy) 
MSA ID 1 
z=11.416 


[] 
i| 
[] 


(CEERS-93316) 
MSA ID 28 MSA ID 355 MSA ID 0 
z = 8.763 z = 7.925 z = 4.912 


(CEERS-DSFG-1) i 
MSA ID 2 L3 MSA ID 2763 MSA ID 3149 
z —4.910 z = 4.902 + + z = 4.901 


Extended Data Figure 3: Montage of images of the galaxies targeted for 
MSA spectroscopy presented here. Each frame shows a 3” x 3” stamp centered 
on the location of each galaxy, with ID labeled. The red-green-blue image for 
each galaxy corresponds to the JWST /NIRCam F444W, F356W, and F277W, 
respectively, rotated with North up and East to the left. T'he rectangles show 
the approximate location of the JWST/NIRSpec MSA slit positions. 


and 355. For MSA IDs 1 and 10 we also detect single emission lines that we 
identify as the unresolved [O r1] doublet at redshifts that are consistent with 
(but more precise than) the Lya break redshifts, and we adopt those values as 
our fiducial redshifts for those galaxies. Source ID 69 shows a weak break at 
~1.3 um, and this methodology does find a maximum likelihood of z = 9.744. 
However, the result has a large uncertainty, with a significant probability of 
the redshift being lower. We thus do not consider MSA ID 69 confirmed. We 
present our redshift measurements in Extended Data Table 1 and the spectra 
of the confirmed high-z objects are shown in Extended Data Figure 2. 
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Extended Data Figure 4: 2D and 1D NIRSpec spectra of CEERS-93316, 
CEERS-DSFG-1 and two other galaxies at similar redshifts z ~ 4.9, highlight- 
ing a potential overdensity in the field at this redshift. The two horizontal 
yellow lines in the 2D spectra indicate the extraction aperture used to measure 
the 1D spectra. Dotted purple vertical lines indicate the location of several 
emission lines when present. The shaded grey area corresponds to lo uncer- 
tainties in the 1D spectra and the zero flux level is marked by a dashed black 
horizontal line. 
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For Maisie's Galaxy, CEERS2.588, and MSA ID 64, we considered the 
alternative interpretation of a Balmer break at redshifts zgg = 3.14, 3.02 and 
2.70, respectively. Averaging the observed spectrum in f, units over wave- 
length intervals (1 + zgg)x 3145-3563 A and 3751-4198 A [14], for Maisie’s 
Galaxy we measure S/N of 0.0 and 9.7 below and above the break, and a 30 
limit on the flux ratio > 2.6, which is considerably larger than the maximum 
value of 1.8 that would be predicted from stellar population synthesis mod- 
els [14]. This supports the higher redshift interpretation of the break resulting 
from Lya opacity of the IGM. For CEERS2_588 the spectrum is negative over 
the bluer window; the background may be slightly oversubtracted due to the 
presence of a large galaxy near the offset shutter in this two-shutter slitlet 
(see Extended Data Figure 3). The 3c limits on the break amplitudes for 
CEERS2.588 and 64 are > 2.2 and > 1.9, respectively, also inconsistent with 
the Balmer break interpretation. 

The derived spectroscopic redshifts are in broad agreement with the pho- 
tometric redshift. For Maisie's Galaxy, in particular, the photometric redshift 
has been modified since the galaxy was discovered as data processing improved. 
The galaxy was originally thought to be at z ~ 14, but the photometric redshift 
was revised to z = 11.8*02 after improved astrometric registration recovered 
weak flux in F150W [5]. Improved NIRCam zero-point flux calibration revised 
this further to z = 11.5*02 [9], while the latest photometric redshift from the 
photometry described in this paper (which is updated relative to earlier work 
[9] to better match HST and JWST photometry) suggests z = 11.08*2 75, 
consistent within the errors with the spectroscopic redshift z — 11.416. For 
CEERS2_588 the photometric (z = 11.02* 032) and spectroscopic (z = 11.043) 
redshifts also agree very well. 

The Lya break redshift we measure for CEERS2.588 (z = 11.43*0 01) 
is significantly larger than the [Om] emission line redshift. A smaller offset 
with the same sense is found for Maisie's Galaxy. The spectral breaks in both 
galaxies are more gradual or “rounded” than is typical for Lya breaks observed 
in galaxies at lower redshifts, which may indicate the presence of a strong Lya 
damping wing due to high column densities of neutral gas either intrinsic to (or 
immediately surrounding) the galaxies [47] or from a highly neutral IGM [48]. 

CEERS-93316 and the millimeter-detected galaxy CEERS-DSFG-1 (MSA 
ID 2) both have redshifts z = 4.91. We identify two more galaxies (MSA 
IDs 2763 and 3149) at nearly the same redshift in the NIRSpec observations 
(Extended Data Figure 4, Extended Data Table 1), supporting the hypoth- 
esis of a physical overdensity [16]. Moreover, the NOEMA-detected galaxy 
1.5 arcsec from CEERS-93316 (Extended Data Figure 1) has a very simi- 
lar photometric redshift and may well also be associated, but still requires 
spectroscopy. 


Emission lines 


The spectra of Maisie’s Galaxy and CEERS2_588 exhibit [O 11] emission lines 
with S/N = 4.9 and 6.4, respectively. While not highly significant, the emission 
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lines are spatially well-aligned with the galaxies’ continuum emission, and care- 
ful inspection of subsets of data from the nodded exposures show no indication 
that these features are data artifacts [21]. At the expected wavelengths of the 
[Ne 1113870 line, CEERS2.588 shows a feature with S/N = 2.0, and Maisie’s 
Galaxy shows positive signal with still lower significance. Weak positive fea- 
tures are seen near the expected wavelengths for O 11111660 and C 111|1909, but 
these cannot be considered significant. 

'The emission lines in CEERS-93316 were fitted manually from the spec- 
trum by first subtracting a local continuum and fitting Gaussian profiles to 
the visible emission lines. In the case of the blended but partially-resolved 
[O 111|4960,5008 lines, a double Gaussian function was fitted. Fluxes and upper 
limits for the main emission lines present for the sources in this work are pre- 
sented in Extended Data Table 2. The Balmer lines (Ha and Hf) are corrected 
for underlying stellar absorption using the best-fitting stellar continuum from 
the Bagpipes fits (corrections « 196; see below for a description of the Bagpipes 
fits). 

Ha is blended with [N 11]6550,6585 in the NIRSpec data. We apply a correc- 
tion [N r1]/Ha ~ 2% based on an Oxygen abundance of 12 + log;,(O/H) = 8.0 
estimated from the [O 11], [O11], H8 and [Ne iri] lines. The Ha flux, corrected 
for [Nt] emission and for slit losses, is (5.27 + 0.11) x 10718 erg cem? s71, 
and the Balmer decrement line flux ratio is f(Ho)/ f (H8) = 5.2 -£ 0.5. Assum- 
ing Case B recombination at T = 10* K and Galactic extinction [49] with 
Ry = 3.1, this corresponds to nebular extinction E(B — V) = 0.60 + 0.10 
and a factor of 4.1* 5$ attenuation at Ha. Assuming the nebular emission 
is produced by young star formation, this implies an extinction-corrected 
SFR(Ho) = 28*5 Mo yr-! [50] for a Chabrier IMF [51]. 


Spectral energy distribution modeling 


We use three spectral energy distribution (SED) fitting codes to derive physical 
properties for Maisie's Galaxy, CEERS2.588, and CEERS-93316. SED fitting 
results are summarised in Extended Data Table 3. 

For the two galaxies with z > 11 the three methods yield results for the 
stellar masses that are consistent within the uncertainties, and within a factor 
of about two for the best-fitting values. Dust attenuation values (Ay) are 
also consistent within the uncertainties, although Cigale finds a larger best- 
fit value with larger uncertainties. The SFRs span larger ranges (factors of 3 
to 4 between methods), likely due to different assumptions on allowable star 
formation histories (SFHs) and, for Cigale, the larger derived value for dust 
attenuation. 

For CEERS-93316 the parameters derived by the three methods span larger 
ranges, most likely reflecting the complexity of this galaxy with its unusual 
combination of strong nebular line emission and large dust reddening. Bag- 
pipes finds the lowest values for M, and Ay, while Cigale and Synthesizer 
find higher values that are consistent within their uncertainties. The large 
reddening implied for both the stars and the ionized gas is atypical for low 
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Extended Data Table 3: Summary of galaxy physical properties 


MSA ID Stellar Mass SFR Stellar Ay 
logio(Ms/Mo) (Mo yr~") (mag) 

Bagpipes [54] 1 8.4403 32 0,070.05 
10 8.5704 42 ooze tt 

0 8.4+0.1 341 1.2 + 0.4 

Cigale [60, 61] 1 8.7 + 0.3 8 0.4 + 0.4 
10 8.8 + 0.3 1205 0.4 + 0.3 

0 9.0 + 0.2 18+ 11 2.1 + 0.6 

Synthesizer [66, 67] 1 8.6 + 0.3 22 1 0.1 + 0.1 
10 8.7 + 0.1 10+ 4 0.2 + 0.1 

0 8.8 + 0.1 60 + 20 2.3 + 0.2 


Extended Data Table 3: Stellar population properties for Maisie’s Galaxy 
(MSA ID 1, z = 11.416), CEERS2_588 (MSA ID 10, z = 11.043) and CEERS- 
93316 (MSA ID 0, z = 4.912) from the three different codes described in the 
Methods section. All uncertainties listed are 68% confidence intervals. 


mass (M, < 10? Mc) galaxies at high redshift [52, 53]. Bagpipes and Synthe- 
sizer find the lowest and highest values for the SFR, respectively, differing by 
a factor of 20, while Cigale’s best fit SFR is intermediate between the other 
methods with a larger fractional uncertainty. The larger reddening and SFR 
values from Cigale and Synthesizer may arise because their fits are addition- 
ally constrained by the emission line EWs (and, with Synthesizer, the spectral 
continuum slope), whereas line emission influences Bagpipes only through its 
contribution to the NIRCam photometry. Bagpipes also permits greater flexi- 
bility in the wavelength dependence of dust attenuation than was allowed for 
the other methods. 

The SFR derived from the extinction-corrected Ha emission from CEERS- 
93316 (287% Mo yr-!) is on the high end of the range derived by Cigale, and 
below the range derived by Synthesizer. The evidence for strong reddening 
suggests the possibility that some star formation takes place in regions that 
are optically thick at UV-optical rest frame wavelengths, and hence that SFRs 
from SED fitting and/or Ha emission may be underestimated. The 30 upper 
limit at 1.1 mm from NOEMA sets an upper limit to the dust-obscured SFR 
of ~ 20-200 M yr-! assuming modified black body emission with dust tem- 
peratures in the range 74 = 40-70 K and a typical spectral index of a = 1.8. 
These limits are reasonably consistent with the SFR estimates derived from 
SED modeling and Ha, although the Ha and Synthesizer SFRs would imply 
dust temperatures Ty Z 41 K and = 48 K, respectively. 

'This analysis assumes that the light from CEERS-93316 is dominated by 
emission from stars and from gas ionized by star formation. Active galac- 
tic nuclei (AGN) can also be responsible for strong line emission, and indeed 
JWST observations have identified a significant population of reddened, broad 
line AGN at z ~ 5 [27, 28]. Most of these AGN are very compact or unresolved, 
whereas CEERS-93316 is spatially extended in NIRCam imaging, indicating 
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that galaxy starlight is present, but this does not exclude a possible contribu- 
tion from AGN emission as well. There is a tentative emission feature at the 
wavelength of Mg 11]2796,2804, a doublet frequently detected in AGN, but its 
significance is marginal (S/N — 1.8). 


Bagpipes 


The Bagpipes [54] fits make use of stellar population synthesis models [55] 
assuming a Kroupa IMF [56] with mass range 0.1-100 Mo, which results in 
masses and SFRs very similar to those from the Chabrier IMF [51] used for 
the Cigale and Synthesizer fits. We consider stellar metallicities in the range 
0.01 € Z/Zo < 3.5. Gas metallicity is set equal to that of the stars. Nebular 
emission lines and continuum are included in the model using results from 
the Cloudy photoionization code [57]. Dust attenuation is included using a 
flexible-slope model, which also includes a Drude profile to model the 2175 A 
dust bump and an additional factor multiplying the attenuation applied to 
star-forming regions (which we take to include all stellar emission arising from 
stars younger than 10 Myr and all resulting nebular emission) [58]. The details 
of further Bagpipes model assumptions are identical to those used in previous 
work [6, 59]. In particular, we use a constant SFR model, shown to produce 
stellar ages and stellar masses at the lower end of the plausible range for 
objects of interest. We apply a Gaussian prior to the redshift of each object 
based on the new spectroscopic data, with a standard deviation of 0.01 in each 
case. The median value for the Gaussian prior is set to z = 11.42 for Maisie’s 
Galaxy, z = 11.04 for CEERS2_588, and z = 4.91 for CEERS-93316. For the 
Bagpipes fits, we use only galaxy photometry and spectroscopic redshifts, and 
do not include additional constraints from the spectra, such as emission line 
strengths. 


Cigale 


We use the latest version of Cigale [60, 61], using stellar population synthe- 
sis models [62] assuming a Chabrier IMF [51] with mass range 0.1-100 Mo. 
Stellar metallicities are allowed to vary from 0.0007 € Z/Zo < 1 for the 
z > 11 galaxies and 0.007 € Z/Zo < 1 for CEERS-93316. Gas metallicites 
are allowed to vary independently from those of the stars, with a minimum 
value 10x larger than the stellar minimum metallicity. We use updated emis- 
sion line models computed from Cloudy [57]. A Calzetti extinction law [63] is 
adopted for the dust attenuation of the stellar continuum. However, the neb- 
ular emission (continuum and lines) is attenuated with a screen model and a 
Small Magellanic Cloud (SMC) extinction curve [64]. No 2175 A dust bump 
is added. For the z > 11 galaxies we that nebular extinction exceeds stellar 
extinction by a constant factor [63], while for CEERS-93316 we assume that 
young stars (< 10 Myr) and ionized gas are more heavily attenuated than older 
stars with a variable scaling factor [65]. We select a “delayed + burst” SFH. 
This form includes a delayed exponential term with SFR œ t xexp(—t/7), with 
T = 200 Myr for the z > 11 galaxies and 500 and 1000 Myr for CEERS-93316. 
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On top of this, a short burst is allowed just prior to the time of observa- 
tion. For the z > 11 galaxies the burst duration is 1 Myr and the burst 
mass fractions considered are 0%, 0.5%, and 10%, whereas for CEERS-93316 
we consider burst durations 5, 10 and 20 Myr and burst mass fractions 5%, 
7.5%, 10%, 12.5% and 15%. For CEERS-93316, the best-fitting mass fraction 
is 10+3.5%; the bursts are smaller for the z > 11 galaxies. The SFR reported 
here is the value averaged over 10 Myr prior to the time of observation. The 
values for Gaussian priors on redshift were the same as for the Bagpipes fit- 
ting. Fitting with Cigale is constrained by the galaxy photometry and redshift. 
Additionally, for CEERS-93316, we use constraints from the EWs measured 
from the NIRSpec spectrum for Ho, Hf, and the two [O i11] lines. The EWs 
from the best-fitting models match the observed values within 3096. The EW 
of the blended [O ri] doublet was included when fitting Maisie's Galaxy and 
CEERS2.588. 


Synthesizer 


With Synthesizer [66, 67], we use stellar population synthesis models [55] 
assuming a Chabrier IMF [51] with mass range 0.1-100 Mc. We sample stellar 
metallicities in the range 0.005 < Z/Zo < 1. Nebular continuum and hydro- 
gen and helium emission lines are included in the modeling, assuming Case B 
recombination, gas temperature 104 K, density 10? cm-?, Lyman continuum 
photon escape fraction 1596, Lya forest transmission 1096 [67]. Conventional 
assumptions about the attenuation law and differential nebular-stellar dust 
effects are assumed [63], implying that the stellar continuum is less reddened 
than the nebular emission by a factor of ~2. Attenuations are allowed to vary 
between Ay = 0 and 1 mag for the z > 11 galaxies, and up to 5 mag for 
CEERS-93316. We explore star formation histories described by a delayed 
exponential form (as also used for Cigale) with timescales ranging from 1 Myr 
to 100 Myr in 0.1 dex steps. The SFR reported here is the value averaged 
over 100 Myr prior to the time of observation. Redshifts are fixed to the val- 
ues given in Extended Data Table 1. Fitting with Synthesizer is constrained 
using both the galaxy photometry and the spectral fluxes, as well as the mea- 
sured redshift. Additionally, for CEERS-93316, Ha and HG EWs measured 
from the spectra are included as a fitting constraint, but all models return 
lower values than those measured by a factor of ~ 4. For Maisie’s Galaxy, a 
relatively strong degeneracy in the age is observed depending on the strength 
of the Balmer break, which is constrained by the spectroscopic data with rel- 
atively low S/N. Fitting spectral data points bluewards of the Balmer break 
yields younger mass-weighted ages and smaller stellar masses (30 + 10 Myr, 
log; (M. /Mc) = 8.30.1). Including NIRSpec spectroscopic data points red- 
wards of the Balmer break increases the derived ages and masses (1204-30 Myr, 
log; (M. /Mc) = 8.8 + 0.1). 
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Photometric redshift mimicry 


The combination of very high EW nebular line emission, strongly reddened 
stellar continuum, and a very specific redshift value is responsible for CEERS- 
93316 mimicking the colors of a galaxy at z ~ 16. At 4.88 S z < 5.09, Ha 
falls in the “triple overlap zone” for all three of the reddest NIRCam band- 
passes used in CEERS (F444W, F410M and F356W), while [O 11|+H/ lies 
within the next bluest filter (F277W), as shown in Figure 3. The emission 
lines boost the signal in those four bands to emulate a roughly flat f, color. 
The highly reddened galaxy continuum fades at rest frame UV wavelengths 
(Aops < 2.4 um), without significant contribution from other emission lines, 
leading to a faint detection in F200W and no significant detections (at CEERS 
exposure times) in the bluer NIRCam filters or by HST. This leads to the 
mistaken impression of a very distant Lya break “dropout” (Figure 4). Galax- 
ies with very strong line emission are common at high redshift, but here the 
combination of strong line emission with a very red continuum is key to fool- 
ing the photometric redshift codes. A similar phenomenon was recognized for 
the submillimeter galaxy CEERS-DSFG-1 (MSA ID 2) [18], which has possi- 
ble (if not necessarily favored) photometric redshift solutions at zphot > 12, 
although the weak [O 111] emission in that galaxy softens the effect compared 
to that for CEERS-93316. The other two galaxies with NIRSpec-confirmed 
redshifts z ~ 4.9 (MSA IDs 2763 and 3149) have bluer colors and are well- 
detected at NIRCam F115W and HST F814W, and therefore do not have 
viable photometric redshift solutions at zynoc > 6. Galaxies like CEERS-93316 
at other redshifts could produce similar photometric aliasing effects, depend- 
ing on where their emission lines fall with respect to photometric bandpasses. 
For example, at 6.71 < z < 6.99 the [O 1115008 line would fall in the same NIR- 
Cam triple-overlap zone for F356W+F410M+F444W, mimicking the colors 
expected at z = 20. Observations through additional intermediate band filters 
could help to further distinguish low and high photometric redshift solutions 
for galaxies like these, but spectroscopy is the ultimate arbiter. 


Data Availability 


The JWST NIRSpec data are available from the Mikulski Archive for Space 
Telescopes (MAST; http://archive.stsci.edu), under program ID 2750. The 
CEERS JWST imaging data are available from MAST under program ID 
1345. Reduced NIRCam data products from the CEERS team are available at 
https: / /ceers.github.io. The NOEMA data is part of program D22AC, which 
will be available at http://vizier.cds.unistra.fr/viz- bin/ VizieR-3?-source-B/ 
iram/noema. 


Code Availability 


JWST NIRSpec data were reduced using the JWST Pipeline (Version 1.8.5, 
reference mapping 1069; https://github.com/spacetelescope/jwst). NIRSpec 
data inspection used the Mosviz visualization tool (https://jdaviz.readthedocs. 
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io/en/latest/mosviz/index.html [43]). Photometric redshifts were measured 
using EAZY [39] (https: //github.com/gbrammer/eazy-photoz/). Emission line 
profiles were fitted and line fluxes measured using LIME [68] (https:// 
lime-stable.readthedocs.io). Three different codes were used for stellar pop- 
ulaation synthesis modeling of galaxy spectro-photometric data: Synthesizer 
[66, 67]; Bagpipes [54] (https: //bagpipes.readthedocs.io); and Cigale [60, 61] 
(https: //cigale.lam.fr). For this paper, we use a non-public version of Cigale 
that fits photometric and spectroscopic data together. 
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